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Direct observation of non-classical photon statistics in parametric downconversion 
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We employ a high quantum efficiency photon number counter to determine the photon number 
distribution of the output field from a parametric downconverter. The raw photocount data directly 
demonstrates that the source is nonclassical by forty standard deviations, and correcting for the 
quantum efficiency yields a direct observation of oscillations in the photon number distribution. 
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Quantum optics 0, Q] is concerned with optical phe- 
nomena that cannot be described by a classical field 
treatment. Operationally, photon correlations are mea- 
sured and nonclassicality is established by violations of 
correlation inequalities imposed by the assumption that 
the field is classical. Well known examples are obser- 
vations of sub-Poissonian photon number statistics 
photon antibunching 0, and photocurrent fluctuations 
below the shot noise level for squeezed light sources 0. 
Although an equivalence between sub-Poissonian statis- 
tics and photon antibunching has sometimes been as- 
sumed, the two phenomena are distinct revelations of a 
nonclassical field 0, , and the choice of test for nonclas- 
sicality may need to depend on the source. Both mea- 
sures are adequate to test nonclassicality for resonance 
fluorescence, but antibunching is certainly not sufficient 
to identify nonclassicality for two-photon sources such 
as an atomic cascade or parametric downconversion 
(PDC) Q. Sophisticated demonstrations of nonclassi- 
cality for two-photon sources have been employed such 
as showing complementarity |lfj| or violating a Bell in- 
equality |ll| . However, direct measurement of the pho- 
ton number distribution may be sufficient to demonstrate 
nonclassicality, and such a measurement is appealing be- 
cause the nonclassicality criterion is strictly based on the 
actual measurement record, inequalities that arise from 
a classical field treatment, and no assumptions whatso- 
ever about the source, propagation or pre-photodetection 
processing. 

Direct measurements of the photon number distribu- 
tion have been elusive for technical reasons, but pho- 
ton counting technology, namely the visible light photon 
counter (VLPC) E|lJ, now provides a means for di- 
rectly testing nonclassicality for the (bunched) output 
field from a PDC as we demonstrate here. Although cri- 
teria for proving nonclassicality of the source via direct 
counting have been suggested, such as sub-Poissonian 
statistics 0|, observing photon number oscillations Q 
and deciding if these oscillations are classical or nonclas- 
sical |ljj, or by employing Hillcry's two criteria of com- 
paring the total probabilities of even vs odd photon num- 
bers being detected 0] , such criteria are not necessarily 



practical, especially for a super-Poissonian two-photon 
source. We introduce a new and practical criterion for 
evaluating nonclassicality for weak bunched sources of 
light and show that our PDC source strongly violates 
this criterion, whereas it does not violate Hillery's in- 
equalities. Our violation of nonclassicality is the first 
time that nonclassical light has been unambigously ob- 
served by direct photon number detection without any 
need for assumptions about the source. 

Furthermore the relatively innocuous assumption that 
detector inefficiency is due to linear loss processes enables 
us to reconstruct from raw data a photon number dis- 
tribution that exhibits even-odd photon number oscilla- 
tions, which strongly violate Hillery's criteria. Although 
photon number oscillations have been reported these 
have in fact been inferred by performing optical homo- 
dyne tomography and thus involve many assumptions to 
process the data; in contrast our reconstruction of pho- 
ton number oscillations requires no assumptions about 
the source and only one assumption about the loss mech- 
anism for detectors. 



Test for non-classical statistics 

With an ideal (i.e., perfect efficiency) photon number 
detector, direct observation of PDC output is predicted 
to exhibit photon number oscillations 2]. For example, 
with P n the probability of observing n photons, an ideal 
detector should yield P\ — P3 = and P2 > 0. How- 
ever, an imperfect detector described by a linear loss 
model may not directly observe photon number oscilla- 
tions, and thus we must establish a criteria for photon 
number statistics to be non-classical. Define T as the 
ratio 



Pi 



P1+P2 + P3 



(1) 



For a perfect detector measuring PDC output, we predict 
r = 1 ; detector inefficiencies will lead to a reduced value 
of this ratio. 
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We now prove that any semi-classical theory of light 
which is constrained to distributions of Poissonian pho- 
ton statistics cannot yield states with T greater than a 
maximum classical bound r c iassicai ■ For a Poisson photon 
number distribution given by P n {n) = c - ™^, this ratio 
has a maximum value r c i ass i ca i = 3/(3 + 2^6) — 0.379, 
saturated by the Poisson distribution with average pho- 
ton number h — y/6. However, one can show that this 
optimal value holds not only for a Poisson distribution, 
but for any weighted sum of Poisson distributions. Con- 
sider a weighted sum P n = aP™ ax + (1 — a)P^ of two 
Poisson distributions P™ ax and P' n , where P™ ax has av- 
erage photon number n = \^6, and P' n is any other Pois- 
son distribution. The ratio T for this weighted sum is 



BBO if 

266nm 



2m A/2 

lens plate PBS Atten 




signal 



prism 

Boxcar 
integrator 



Nd:YAG 
4 th Harmonic 
266nm 



r = 



a pmax 



(1 - a)P^ 



a{P[ 



tiinx _|_ pmax _|_ ^Drnax j 



(l-a)(P{ + P^ + P0 
(2) 

Because P™ ax maximizes T for any single Poisson distri- 
bution, the mathematical relation 

x' x ax + (1 — a)x' x 
¥ V ay+(l- a)y' y 

proves that T < 3/(3 + 2-n/6)- Thus, no sum of Poisson 
distributions can give rise to a distribution with T > 

^classical 

All classical light fields will lead to statistics that can 
be expressed as weighted sums of Poisson photon number 
states. Thus, the classical theory of light predicts that 
the inequality 



r < 



3 + 3V6' 



(4) 



cannot be violated. In contrast, one expects that light 
from PDC will lead to a violation of this condition, which 
can be demonstrated by simply measuring Pi, Pi, and 
Ps- 

In the presence of imperfect detection efficiency, loss 
may serve to degrade this ratio and a violation of the 
classical criterion may not be observed. Consider a PDC 
experiment in which the pump is sufficiently weak that 
the probability of generating more than one photon pair 
is very small. In this case the ratio in Eq. is given 
by r = 77/(2 — 77), where r\ is the detection efficiency. A 
violation of the inequality is not predicted unless 77 > 
3/(3+ VE) « 0.55. 

Fortunately, the Visible Light Photon Counter 
(VLPC) has the capability to detect photon number 
states with high quantum efficiency [H|l3. The VLPC 
has been shown to have quantum efficiencies approaching 
90%. Furthermore, if more than one photon is incident 
on the detector surface, the height of the output electrical 
pulse is proportional to the number of incident photons. 
This gives us information about the number of photons 
that have been detected. 



FIG. 1: Schematic of experimental configuration. The 
BBO crystal is tilted for Type I collinear degenerate phase- 
matching. 



Observation of non-classical statistics 

The experimental setup is shown in Fig. ^ We use 
the fourth harmonic (266nm) of a Q-Switched ND:YAG 
laser as a pump source. The pulse duration of the laser 
is 20ns. Using a pulsed pump allows us to eliminate the 
detector dark counts (20, 000s _1 ) by temporal gating. 
The laser pumps a BBO crystal set for collinear degener- 
ate Type I phase matching (optic axis 47.6 degrees from 
the pump). In this configuration, the down-converted 
photons have half the energy of the pump (532nm), and 
travel in the same direction. The pump is removed by a 
prism, while the down-conversion is focused by a 250mm 
lens onto the detector surface. We also have the option of 
directly illuminating the detector with second harmonic 
light (532nm) from the laser, which is a classical light 
source. The detector output is amplified, and then sent 
to a gated boxcar integrator, which is triggered by the 
laser. The boxcar integrates the pulse over a 20ns win- 
dow, and the output is sent to an Analog-to-Digital con- 
verter and stored on a computer. 

The output of the detector illuminated by light from 
the second harmonic of the laser is shown in Fig. |5J The 
pulse area spectrum features a series of peaks represent- 
ing the different photon number state detections. In 
the inset we show the probability distribution, which is 
calculated by fitting each peak to a Gaussian function. 
The area under each Gaussian curve gives the number 
of events representing that photon number. The area of 
each peak can be normalized by the total area to give the 
probability distribution, which is a Poisson distribution 
as expected. 

Fig.^a) shows the pulse area histogram when the de- 
tector is illuminated by parametric down-conversion, us- 
ing a pump power of \[iW . At this weak pump intensity, 
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FIG. 2: Pulse area histogram from the boxcar integrator ob- 
tained by illuminating attenuated second harmonic (532nm) 
light from the laser onto the detector. The photon number 
distribution (shown in the inset) is a Poisson distribution. 
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FIG. 3: Observation of non-classical statistics, (a) The pulse 
area histogram for detection from the VLPC illuminated by 
parametric down-conversion. The photon number distribu- 
tion is shown in the inset. From this distribution F — 0.442, 
which violates the classical limit by 40 standard deviations, 
(b) Plot of P as a function of pump power. At low pumps 
P drops due to dark counts. At high powers it drops due to 
four photon events, which enhance the three photon proba- 
bility after detection losses. 



a single pump pulse will usually generate zero photons, 
while a photon pair is generated with a small probabil- 
ity. The probability of generating more than one photon 
pair is very small. The figure focusses on the 1, 2, and 
3 photon detection peaks, which we will use to verify 
non-classical statistics. The photon number probability 
distribution is calculated by fitting each peak to a Gaus- 
sian function. These areas are normalized by the total 
area of all the peaks. The calculated probability dis- 
tribution is shown in the inset. One can see that the 
probability of 1 and 2 photon detection is nearly equal, 
but the probability of 3 photon detection is nearly zero. 
These probabilities are Pi = 0.0818, P 2 = 0.0696, and 
P 3 = 0.0061, which yields T = 0.442, representing a 40 
standard deviation violation of the classical limit. 

The large 1 photon probability is due to losses from 
the detector and collection optics. In the limit of low 
excitation, the 1 photon and 2 photon probability can be 
used to calculate the detection efficiency, given by 



From the measurements, it is calculated that the detec- 
tion efficiency is 0.67. Using the measured VLPC quan- 
tum efficiency of 0.85, the photon collection efficiency is 
calculated to be 0.79. 

Fig. OJb) shows the measured value of T as a function 
of pumping intensity. The black line represents the clas- 
sical limit, which is violated for a large range of pumping 
intensities. At high pumping intensities T begins to drop. 
This drop is due to an increase in the two pair creation 
probability, which, in the presence of losses, will enhance 
the 3 photon detection probability. The parameter T also 
drops at low pumping intensities. This drop is attributed 
to the dark counts of the VLPC. At low pumping inten- 
sities the relative fraction of detection events originating 
from dark counts becomes large. This enhances the 1 
photon probability, which reduces the value of V. 

Reconstruction of even-odd oscillations 

With high detector efficiency, the emitted output of 
PDC is predicted to feature even-odd photon number 
oscillations due to the two photon nature of the pro- 
cess. These oscillations lead to the non-classical statistics 
discussed in the previous section. Direct observation of 
these oscillations using the photon counting capability of 
the VLPC would be a remarkable achievement; unfor- 
tunately, direct observation of these oscillations requires 
extremely high quantum efficiencies. Even the relatively 
high detection efficiencies of 0.67 in our experiment are 
not predicted to observe this oscillatory behavior. The 
requirement of very high quantum efficiency makes direct 
observation of the photon number oscillation extremely 
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FIG. 4: Pulse area histograms in (a), (b), and (c) are shown for 4, 6, and 8/j.W pumping powers respectively. The probability 
distributions are shown in the insets. By correcting for the detection efficiency and dark counts we reconstruct the original 
photon number distributions shown on the right. These distribution show the clear even-odd photon number oscillations from 
parametric down- conversion. 



difficult in practice. However, one can make an accu- 
rate independent measurement of the photon detection 
efficiency, and correct for this effect in the photon num- 
ber distribution. This allows the reconstruction of the 
original even-odd oscillations of the field. 

The detection efficiency can be corrected for as fol- 
lows. Define pi as the probability that the photon field 
contained i photons, and fi as the probability that i pho- 
tons are detected. In the presence of losses, these two 
distributions are related by 

fi X (')''' M '<> ><• ( 6 ) 

j=i v ' 

In order to calculate pi from /, , the above transformation 
must be inverted. To perform this inversion, we truncate 
the photon number distribution at some photon num- 
ber n, which is sufficiently large such that p n +i ~ is 
a good approximation. Under this approximation, the 
initial and final probability distributions are simply re- 
lated by a matrix, whose coefficients are given by Eq. JSJl . 



The dark counts of the VLPC can also be accounted for 
by this matrix, which can then be inverted to calculate 
the photon number distribution of the field. It is impor- 
tant to emphasize that there are no fitting parameters 
in this model. The only two parameters, the quantum 
efficiency and dark counts of the VLPC, are both in- 
dependently measured. Once they are known there is a 
one-to-one relationship between the actual and measured 
photon number distribution. 

Fig. 01 shows the result of the photon number recon- 
struction. Three different pumping intensities are shown. 
For each pump intensity, the left panel shows the pulse 
area histogram, and the inset to the panel shows the cal- 
culated photon probability distribution. The right panel 
shows the reconstructed photon number distribution us- 
ing the measured quantum efficiency and backgrounds. 
The photon number distribution is truncated at 10 pho- 
tons. The reconstructed probabilities demonstrate very 
clear even-odd oscillations as predicted from PDC. 

At higher photon numbers, it can be seen that the re- 
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constructed distribution becomes slightly negative. This 
erroneous effect is caused by truncation error. As the 
pumping intensity is increased, the approximation that 
the photon distribution can be truncated after 10 pho- 
tons becomes less accurate. This error manifests itself 
in the probabilities becoming slightly negative for the 
9 and 7 photon probability. This error is worst at the 
largest pumping intensity of 8/jW, where the truncation 
approximation is least accurate. One could suppress this 
error by truncating at a higher photon number. Unfor- 
tunately, because of the limited range of the amplifiers 
and A2D converters, it is difficult to measure these higher 
order photon numbers in practice. This puts a limit on 
the pumping power one can use and still get a good re- 
construction. It is possible that an improved numerical 
algorithm over simply putting a cutoff in the number 
distribution may overcome some of these practical diffi- 
culties. 



Conclusion 

In conclusion, we have directly observed non-classical 
photon counting statistics from PDC. We have shown 
theoretically that the photon counting statistics for all 
classical fields must satisfy the inequality given in Eq. 
Using the high quantum efficiency and photon number 
detection capability of the VLPC we have experimen- 
tally demonstrated violations of this inequality by light 
emitted from PDC. By correcting for the quantum ef- 



ficiency and dark counts of the VLPC, we have also 
succeeded in reconstruction the even-odd oscillation in 
the photon number distribution of light generated by the 
down-conversion field. 
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